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ABSTRACT 

We consider stochastic reconnection in a magnetized, partially ionized 
medium. Stochastic reconnection is a generic effect, due to field line wander- 
ing, in which the speed of reconnection is determined by the ability of ejected 
plasma to diffuse away from the current sheet along magnetic field lines, rather 
than by the details of current sheet structure. As in earlier work, in which we 
dealt with a fully ionized plasma, we consider the limit of weak stochasticity, so 
that the mean magnetic field energy density is greater than either the turbulent 
kinetic energy density or the energy density associated with the fluctuating com- 
ponent of the field. For specificity, we consider field line stochasticity generated 
through a turbulent cascade, which leads us to consider the effect of neutral drag 
on the turbulent cascade of energy. In a collisionless plasma, neutral particle 
viscosity and ion-neutral drag will damp mid-scale turbulent motions, but the 
power spectrum of the magnetic perturbations extends below the viscous cutoff 
scale. We give a simple physical picture of the magnetic field structure below 
this cutoff, consistent with numerical experiments. We provide arguments for 
the reemergence of the turbulent cascade well below the viscous cut-off scale and 
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derive estimates for field line diffusion on all scales. We note that this explains 
the persistence of a single power law form for the turbulent power spectrum of 
the interstellar medium, from scales of tens of parsecs down to thousands of kilo- 
meters. We find that under typical conditions in the ISM stochastic reconnection 
speeds are reduced by the presence of neutrals, but by no more than an order 
of magnitude. However, neutral drag implies a steep dependence on the Mach 
number of the turbulence. In the dense cores of H2 regions the reconnection 
speed is probably determined by tearing mode instabilities. 

Subject headings: Magnetic fields; Galaxies: magnetic fields, ISM: molecular 
clouds, magnetic fields; Stars: formation 

1. Introduction 

One of the fundamental properties of astrophysical magnetic fields is their ability to 
change topology via reconnection (see Priest & Forbes, 2000) . It is impossible to understand 
the origin and evolution of large scale magnetic fields without understanding the mobility 
of magnetic field lines. In a typical astrophysical plasma, resistivity is very small and fiux 
freezing, which follows from assuming zero resistivity, should be an excellent guide to the 
motion of magnetic fields. In spite of this, fast magnetic dynamo theory (see Parker, 1979; 
Mofi'att, 1978; Krause & Radler, 1980) invokes a constantly changing magnetic field topology 
and motions during which field lines cross each other^. This assumption is supported by 
observations of the solar magnetic field (see Dere, 1996; Innes, Inhester, Axford & Wilhelm, 
1997, and references contained therein) which are difficult to explain unless fiux freezing is 
routinely violated on time scales short compared to resistive time scales, at least within thin 
current sheets. 

The Sweet-Parker model of reconnection is the simplest and most robust (Parker, 1957; 
Sweet, 1958). In this model, reconnection takes place within a thin current sheet, which 
separates two large volumes containing uniform, and very different, magnetic fields. The 
resulting reconnection speed is less than the Alfven speed by the square root of the Lundquist 
number ~ T^l^^"^ — {v/^aLY^^ , where 1] is the resistivity, Va is the Alfven speed, and L is 
the length of the current sheet, assumed to be determined by the large scale geometry of 



^Whilc standard mean-field dynamo theory has been severely criticized on theoretical and numerical 
grounds (Vainshtein & Cattaneo, 1992; Gruzinov & Diamond, 1994, 1996; Cattaneo & Hughes, 1996; Hughes 
et al., 1996; Brandenburg, 2001), an alternative version can be formulated which evades these criticisms 
(Vishniac & Cho, 2001) and which assumes fast reconnection only on two dimensional surfaces. 
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the problem. Under typical astrophysical conditions this is very slow (e.g. for the Galaxy 
as a whole IZl ~ 10^°). This reconnection speed is set by a geometrical constraint. Indeed, 
plasma tied to the reconnecting magnetic field lines must be ejected from the ends of the 
narrow current sheet. The disparity of scales, one of which is macroscopic/astrophysical, 
while the other is microscopic, i.e. determined by ohmic diffusion, makes the reconnection 
slow. 

This evident shortcoming of the Sweet-Parker reconnection scheme has stimulated inter- 
est in alternative models that allow fast reconnection. Although the literature on magnetic 
reconnection is extensive (e.g. Priest & Forbes, 2000, and references therein) it does not 
successfully address this question. Models that invoke an X-point reconnection geometry 
(Petschek, 1964) have been shown to be unstable for sufficiently high (see Biskamp, 
1996), while anomalous resistivity fails to provide rapid reconnection under most astrophys- 
ical conditions (see Parker, 1979). A general review of astrophysical magnetic reconnection 
theory can be found in Bhattacharjee, Ma, & Wang (2003). 

A notable exception to this discouraging state of affairs is the recent work on fast col- 
hsionless reconnection (Biskamp, Schwarz & Drake, 1997; Shay, Drake, Denton & Biskamp, 
1998; Shay & Drake, 1998, see also the discussion by Bhattacharjee, Ma and Wang 2001). 
This work indicates that under some circumstances a kind of standing whistler mode can 
stabilize an X-point reconnection region. However, these studies have not demonstrated the 
possibility of fast reconnection for generic field geometries. They assume that there arc no 
bulk forces acting to produce a large scale current sheet, and that the magnetized regions 
are convex, which minimizes the energy required to spread the field lines. In addition, this 
mechanism requires a coUisionless environment, where the electron mean free path is less 
than the current sheet thickness (Trintchouk, Yamada, Ji, Kulsrud & Carter, 2003; Ji, Ya- 
mada, Hsu & Kulsrud, 1998). Unfortunately, in the laboratory the current sheet thickness 
is comparable to the ion Larmor radius and it is unclear how to generalize this criterion to 
the interstellar medium, where the Sweet-Parker current sheet thickness is typically much 
greater than the Larmor radius. // we require that the ion skin depth, the characteristic 
scale of the standing Whistler mode, be greater than the current sheet thickness, then we 
have a criterion which is rarely satisfied in the interstellar medium. In this paper we will 
concentrate on a mechanism that acts in turbulent media and when it works, produces rapid 
reconnection under a broad range of field geometries, without regard to the particle collision 
rate. We will defer all discussion of the relationship between stochastic reconnection and 
small scale coUisionless effects to a later paper. 

In a previous paper (Lazarian & Vishniac, 1999, henceforth LV99) we discussed 'stochas- 
tic reconnection', a process which is similar to Sweet-Parker reconnection, except that 
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stochastic wandering of the field lines produces a broad outflow region. The properties 
of the outfiow region are insensitive to the width of the current sheet (and the value of TZl), 
but depend on the level of field line stochasticity. In a sufficiently noisy environment the 
reconnection speed becomes a large fraction of the Alfven speed. In an extremely quiet envi- 
ronment the field lines do not enter or leave the current sheet over its entire length and we 
recover the Sweet-Parker reconnection model. In LV99 we dealt with an inviscid and totally 
ionized fluid. ^ 

The notion that magnetic field stochasticity might affect current sheet structures is 
not unprecedented. In earlier work Speiser (1970) showed that in collisionless plasmas the 
electron collision time should be replaced with the time a typical electron is retained in the 
current sheet. Also Jacobson & Moses (1984) proposed that current diffusivity should be 
modified to include diffusion of electrons across the mean field due to small scale stochasticity. 
These effects will usually be small compared to effect of a broad outflow zone containing both 
plasma and ejected shared magnetic flux. Moreover, while both of these effects will affect 
reconnection rates, they are not sufficient to produce reconnection speeds comparable to the 
Alfven speed in most astrophysical environments. 

It is important to distinguish between stochastic reconnection, as discussed in LV99, 
and the more conventional notion of turbulent reconnection. The latter usually involves 
substituting a turbulent diffusivity for the resistivity, and involves a degree of small-scale 
mixing which is forbidden on energetic grounds (see §2 in Parker, 1992, and the references 
contained therein). A common variation of this hypothesis is that instabilities in the current 
sheet will produce a hugely broadened current sheet and a large effective resistivity within it 
(for a review see Biskamp, 2000). On the contrary, the former is largely a topological effect, 
using conventional estimates of resistivity, and the only strong mixing associated with it has 
to do with the polarization of field fines crossing the current sheet. (That is, one expects 
sharp gradients, in the current sheet, in the component of the magnetic lines perpendicular 
to the current sheet.) In this sense, stochastic reconnection belongs to the class of models 
which try to explain fast reconnection by appealing to a current sheet geometry which is 
'natural' in some sense, but evades the limits set by the Sweet-Parker model. (The bulk 
of the discussion in Priest & Forbes (2000) is centered on laminar three dimensional field 
configurations which can lead to similarly rapid reconnection speeds.) 

In a recent paper Kim & Diamond (2001) addressed the problem of stochastic recon- 
nection by calculating the turbulent diffusion rate for magnetic flux inside a current sheet. 



^In this paper we shall show that this approximation is valid for partially ionized collisionless plasma up 
to a certain percentage of neutrals, and for coUisional gases with resistivity larger than viscosity. 
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They obtained similar turbulent diffusion rates for both two dimensional and reduced three 
dimensional MHD. In both cases the presence of turbulence had a negligible effect on the 
flux transport. The authors pointed out that this would prevent the anomalous transport of 
magnetic flux within the current sheet and concluded that both 2D and 3D stochastic recon- 
nection proceed at the Swcct-Parker rate even if individual small scale reconnection events 
happen quickly. If true this would be not only rule out the LV99 reconnection scheme, but 
also any other fast reconnection scheme. In general astrophysical plasmas are turbulent and 
if the enhancement of the local reconnection speed, e.g. due to coUisionless effects (see Drake 
et al. 2001), is irrelevant then reconnection must always be slow. 

However turbulent diffusion rates within the current sheet are irrelevant for the process 
of stochastic reconnection or, for broadly similar reasons, fast coUisionless reconnection. 
The basic claim in LV99 is that realistic magnetic field topologies allow multiple connections 
between the current sheet and the exterior environment, which would persist even if the 
stochastic magnetic field lines were stationary ("frozen in time") before reconnection. This 
leads to global outflow constraints which are weak and do not depend on the properties of the 
current sheet. In particular, the analysis in LV99 assumed that the current sheet thickness 
is determined purely by ohmic dissipation and that turbulent diffusion of the magnetic field 
is negligible inside, and outside, the current sheet. The major uncertainty in this model is 
the behavior of the reconnected flux elements, which are nearly perpendicular to the current 
sheet and must undergo multiple reconnections before being ejected. We note also that 
models of coUisionless reconnection also evade the objection posed by Kim and Diamond 
topologically, that is, by stabilizing an X-point reconnection topology, and opening up the 
rest of the current sheet. 

"Hyper-resistivity" (Strauss, 1985; Bhattacharjee & Hameiri, 1986; Hameiri & Bhat- 
tacharjee, 1987) is a more subtle attempt to derive fast reconnection from turbulence within 
the context of mean-field resistive MHD. The form of the parallel electric field can be derived 
from magnetic helicity conservation. Integrating by parts one obtains a term which looks 
like an effective resistivity proportional to the magnetic helicity current. There are several 
assumptions implicit in this derivation, but the most important problem is that by adopt- 
ing a mean-field approximation one is already assuming some sort of small-scale smearing 
effect, equivalent to fast reconnection. Strauss (1988) partially circumvented this problem 
by examining the effect of tearing mode instabilities within current sheets. However, the 
resulting reconnection speed enhancement is roughly what one would expect based simply 
on the broadening of the current sheets due to internal mixing. This effect does not allow 
us to evade the constraints on the global plasma flow that lead to slow reconnection speeds, 
a point which has been demonstrated numerically (Matthaeus & Lamkin, 1985) and analyt- 
ically (LV99). Nevertheless, we show in §4 that this effect may be important in the densest 
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and coldest parts of the ISM. 

A partially ionized plasma fills a substantial volume within our galaxy and the earlier 
stages of star formation take place in a largely neutral medium. This motivates our study 
of the effect of neutrals on reconnection. The role of ion-neutral collisions is not trivial. On 
one hand, they may truncate the turbulent cascade, reducing the small scale stochasticity 
and decreasing the reconnection speed. On the other hand, the ability of neutrals to diffuse 
perpendicular to magnetic field lines allows for a broader particle outfiow and enhances 
reconnection rates. 

Reconnection in partially ionized gases has been already studied by various authors 
(Naidu, McKenzie & Axford, 1992; Zweibel & Brandenburg, 1997). In a recent study (Vish- 
niac & Lazarian, 1999, henceforth VL99) we studied the diffusion of neutrals away from the 
reconnection zone assuming anti-parallel magnetic field lines (see also Heitsch & Zweibel, 
2003a) The ambipolar reconnection rates obtained in VL99, although large compared with 
the Sweet-Parker model, are insufficient either for fast dynamo models or for the ejection 
of magnetic fiux prior to star formation. In fact, the increase in the reconnection speed 
stemmed entirely from the compression of ions in the current sheet, with the consequent 
enhancement of both recombination and ohmic dissipation. This effect is small unless the 
reconnecting magnetic field lines are almost exactly anti-parallel (Vishniac & Lazarian, 1999; 
Heitsch & Zweibel, 2003b) . Any dynamically significant shared field component will prevent 
noticeable plasma compression in the current sheet, and lead to speeds practically indis- 
tinguishable from the standard Sweet-Parker result. Since generic reconnection regions will 
have a shared field component of the same order as the reversing component, the implication 
is that reconnection and ambipolar diffusion do not change reconnection speed estimates 
significantly. 

None of this work on reconnection in partially ionized plasmas includes the effects of 
stochasticity. We expect that in the presence of turbulence, reconnection rates will be 
substantially enhanced, as they are in completely ionized plasmas. To generalize the concept 
of stochastic reconnection to partially ionized plasmas we need a model for the small scale 
stochasticity of a turbulent magnetic field in a partially ionized plasma. In this paper we 
will begin by considering this problem, and then apply our results to the reconnection speed. 

Following LV99, we consider reconnection in the presence of a weakly stochastic mag- 
netic field. Except for the presence of noise, we imagine a reconnection event exactly like 
a generic Sweet- Parker reconnection event. Two volumes with average magnetic fields that 
are of comparable strength, but differing directions are in contact over a surface of length 
L. Due to the stochastic nature of the fields, field lines come into contact over many small 
patches (see Fig. 1). For each individual patch the Sweet-Parker reconnection model should 
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be applicable (or at least constitute a minimal reconnection speed). The enhancement of 
reconnection rates follows from two effects. First, since individual field lines wander out of 
the narrow current sheet relatively easily, the longitudinal patch size is much smaller than 
the overall size of the system. This reduces the effective value oi TZ^ and raises the local 
reconnection speed. Second, whereas in the Sweet-Parker scheme magnetic field lines recon- 
nect sequentially, in the presence of field line wandering field lines many independent patches 
are brought into direct contact and can reconnect simultaneously. As a consequence, the 
rate of reconnection of the magnetic flux is increased by a large factor, whose exact value 
depends on the level of noise in the system. 

In LV99 the truncation of the turbulent cascade was assumed to be due to resistivity. 
Consequently, the smallest scale of field line wandering decreases as resistivity decreases, and 
the number of independent patches in contact within the reconnection zone increases. Prom 
this we concluded that for an idealized inviscid fluid the reconnection rate does not depend 
on fluid resistivity. It does depend on the level of magnetic field stochasticity. In the specific 
case where the field line stochasticity is caused by a turbulent cascade, it depends on the 
amplitude of the turbulence. However, there is no necessary connection between turbulent 
motions and the reconnection speed. In particular, the rate of turbulent transport of mean 
magnetic fiux is assumed to be negligible in this model, and even the complete absence of 
turbulent diffusion in the current sheet would not reduce the stochastic reconnection rate 
(cf. Kim & Diamond, 2001). 

To quantify stochastic reconnection we have to use a particular description of turbu- 
lence. Motions in a magnetized medium can be expanded into incompressible (Alfven) and 
compressible (fast and slow) modes. There are theoretical arguments (Goldreich & Sridhar, 
1995; Lithwick & Goldreich, 2001) (see Cho, Lazarian & Vishniac 2003a, for a review) sug- 
gesting that the nonlinear cascade of power for these modes proceeds separately, although 
not entirely independently. Simulations in Cho & Lazarian (2002) support this idea and show 
that the Goldreich & Sridhar (1995, henceforth GS95) scaling is applicable to Alfvenic part 
of the MHD cascade. For our purposes the scalings of slow and fast modes (Cho & Lazarian, 
2002), (Cho & Lazarian, 2003) are less important since they are subjected to coUisionless 
damping^ 

Here we assume that the GS95 model describes incompressible turbulence above the 



•^Whether or not those damped modes are important depends on the process studied. For instance, Yan 
& Lazarian (2002) show that for scattering of cosmic rays the residual small amplitude fast modes are much 
more efBcient than the Alfven modes. This, however, is not true for the field wandering that we consider in 
this paper. 
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ambipolar damping scale^. To describe MHD turbulence below the scale of viscous damp- 
ing we present a new model of magnetic field structure in this regime. This model is in 
rough agreement with numerical simulations by Cho, Lazarian & Vishniac (2002b, hence- 
forth CLV02b). 

In §2 of this paper we will consider the effect of a large neutral fraction on a strongly 
turbulent cascade in a magnetized plasma. In §3 we apply this to the problem of reconnection 
in partially ionized plasmas. In §4 we apply this work to various phases in the ISM. Finally, 
§5 contains our basic conclusions. 

2. Magnetohydrodynamic Turbulence in a Pctrtially Ionized Plasma 

In this section we will consider the effect of neutral particles on the turbulent cascade in 
the ISM. We begin by briefly reviewing the nature of the strong turbulent MHD cascade and 
the dynamical influence of neutral particles. In §2.2 we describe the cascade when viscous 
damping, due to neutral particles, is strong, but the one fluid approximation remains valid. 
In §2.3 we consider the uncoupled regime, covering scales where the neutral particles exert 
a uniform drag on all motions. We end, in §2.4, with a brief discussion of the implications 
of this picture for observations of turbulence in the ISM. 

2.1. Neutral- ion damping 

The role of neutral-ion damping in MHD turbulence has been discussed previously in 
the context of the ISM (in particular, see Spangler, 1991; Minter & Spangler, 1997). The 
basic conclusion was that neutral fluid heating is a plausible sink for the turbulent energy 
revealed through measurements of interstellar scintillation. Here we are concerned instead 
with how a neutral gas component will modify the turbulent power spectrum. The most 
relevant observational point is that the ISM turbulent power spectrum has no strong fea- 
tures at wavelengths where neutral- ion coupling would be expected to play a dominant role 
(Armstrong, Rickett & Spangler, 1997). Instead, the power spectrum extends to very small 
scales (< 10^ cm) in an approximate power law. Qualitatively, this suggests that stochas- 
tic reconnection can take place even in partially neutral plasmas. However, several basic 



''We note, however, that our qualitative conclusions for reconnection rates should be valid for other models 
of MHD turbulence (see LV99) as long as they are in rough agreement with observational constraints and 
numerical simulations. 
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questions remain unanswered. Previous work on turbulence in the ISM has not included 
a discussion of the most plausible model for MHD turbulence (although a simple hydrody- 
namic model was addressed, which is remarkably close to the model we use here). Also, we 
need to understand why neutral damping fails to produce a strong signature in the ISM, 
or at least in the diffuse ionized component of the ISM, before we can construct a general 
model for its role in partially ionized plasmas. 



2.1.1. The Goldreich-Sridhar model 

The GS95 model of strong MHD turbulence is based on the notion of a Kolmogorov-like 
cascade with an anisotropy imposed by the large scale magnetic field. The exact degree of 
anisotropy follows from an average balance between hydrodynamic and magnetic forces. Ed- 
dies on a given scale are characterized by a wavenumber perpendicular to the mean magnetic 
field direction, k±, and a parallel wavenumber, k\\, such that the rate of eddy turnover time 
is equal to the rate wave propagation along magnetic field, i.e. 

k±Vk ^ k\\VA , (1) 

where Vk is the typical velocity at the scale characterized by the wavenumber {k\\, k±). As Vk 
goes down with the increase of k± this condition implies eddies that are elongated along the 
field direction, and become more elongated as we go to smaller scales. In order to simplify 
our notation, we will refer to k± below as k. 

If energy is injected isotropically on some scale /, with vi < Va, then the cascade will 
begin in a regime of weak turbulence, in which motions can be characterized as weakly 
interacting waves, with a frequency u = /cyV^ ^constant and a nonlinear decay rate (see 
discussion in a review by Cho, Lazarian & Vishniac 2003) 

so that conservation of energy implies 

Vk ~ vi{kl)-'/\ (3) 

For T~i^ less than the wave frequency the corresponding scale-dependent diffusion coefficient 
will be (LV99) 

Dk ~ ~ ^ (4) 

which is actually independent of scale. The decrease in wave motions at larger k is balanced 
by the decrease in the coherence time of the waves. Here we will ignore the contribution 
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of these scales to turbulent diffusivity in favor of the contribution from smaller, strongly 
turbulent, scales. 

The weak turbulent cascade ends at a scale kr where equation (1) is satisfied. At this 
scale 

vt ~ 77-, and krl ~ ( — 1 . (5) 

Va \vi J 

At larger k the strong turbulence model applies and (see GS95) 

(6) 



The rate of turbulent energy transfer is A;||Va, which means 

while the rms turbulent fluid velocity is given by 

f k\ "^/^ 

v^^vr[^-j . (8) 

The magnetic field perturbations, b^, are 

bk ^ Vk{AT:pf'\ (9) 

This model presupposes that the turbulent velocities are subalfvenic, and we adopt 
this assumption in the rest of this paper. This is less restrictive than it might appear, 
since as long as there is some scale /' in the turbulent cascade where vi' ~ Va we can take 

I ^ rii krp 

, Vt ^ Vl ^ Va and use this model of turbulence for all smaller scales. Moreover, 
if the turbulent energy is larger than the magnetic energy, we can expect rapid growth of the 
magnetic field through the turbulent dynamo, that is, the generation of a disordered field 
in rough equipartition with the turbulent kinetic energy (see, however, Schekochihin et al. 
2003). 

The Goldreich-Sridhar scahngs can be easily understood. They reflect the fact that on 
small scales it is difficult to bend magnetic field lines, but it is still easy to mix them up. 

In a fully ionized astrophysical plasma, shear viscosity is generally less important than 
resistivity in damping MHD turbulence. In a partially neutral medium a combination of 
neutral particle viscosity and ion-neutral coUisional coupling drives damping. Following 
GS95 we concentrate on the diffusion of momentum across field lines. 
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2.1.2. Ion-neutral decoupling: theoretical considerations 

The preceding discussion assumes that the plasma is a single, tightly coupled, fluid 
with negligible viscosity. Many astrophysical fluids are partially ionized, the neutrals are 
imperfectly coupled to the ions. An obvious consequence is a substantially increased viscosity 
since the neutrals can cross magnetic field lines. 

The coupling between ions and neutrals is determined by the rate of ion-neutral colli- 
sions, which is 

t^n = T ""(^^rf^m), (10) 

where is the ion-neutral relative velocity, 0"^ is the ion-neutral collisional cross section, 
rrii and m„ are the typical ion and neutral masses, rin is the neutral number density, and 
angular brackets denote averaging. From (Draine, Roberge & Dalgarno, 1983) we adopt 
{vrCin) ~ 1.9 X 10~^ cm^ s~^. The rate at which neutrals exchange momentum with ions is 

If the mean free path for a neutral particle, /„, in a partially ionized gas with density 
ntot = nn + Tii is much less than the size of the eddies under consideration, i.e. Ink -C 1, the 
damping time is 

-l;.-2 _ . l^tot \ /, \-li-2 



tdamp k ^ I ) i}nCn) k , (11) 

where Vn is the effective viscosity produced by neutrals and c„ is the sound speed in the 
neutrals ^ . 

Consider first a mostly neutral gas (i.e. ^ n„). Turbulence can cascade to small 
scales if the turbulent eddy rate ~ kv}^ is larger than the viscous damping rate i^ampms- 
In a partially ionized gas the one fluid approximation is valid if neutrals experience multiple 
collisions with ions in an eddy turnover time, i.e. > r~^. Therefore, MHD turbulence 
will exhibit Goldreich-Sridhar scaling up to the damping scale if t^/ > t^„^,ping- On the other 
hand, if < t^^^^^-^g neutrals will not follow the ions and magnetic flelds. Instead they will 
form a hydrodynamic cascade as soon as the neutral-ion collisional rate t~l is of the order of 
the eddy turnover time r. At this scale, and all smaller scales, the ionic fluid motions will 



^The viscosity across magnetic field lines, due to ion-ion collisions, is typically small as ion motions are 
constrained by the magnetic field. For a collision rate much smaller than the ion cyclotron frequency the 
ratio of ion perpendicular viscosity to resistivity is a few times the ratio of the ion thermal pressure to the 
magnetic pressure. In the collisional limit this is multiplied by a factor of (Qj/i/j)^ <C 1, and resistivity 
efficiently dissipates magnetic field perturbations on scales greater than the viscous damping scale. The drag 
coefficient for neutral-neutral collisions is ~ 1.5 x 10~^°T^/^ cm^ s~^ with T measured in Kelvins (Spitzer, 
1978), so collisions with other neutrals will dominate for rii/nn less than <^ 0.08T^/^. 
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damp at the rate of ion-neutral collisions , which in mostly neutral gas is much larger 
than t'l. 

What will happen to the magnetic fields when ionic motions are damped? Collisions 
between the ions and the neutral particles will prevent magnetic tension from straightening 
the field lines efficiently. Instead they will be moved, entangled and stretched by large 
undamped eddies. As the eddy turnover rate increases with the decrease of the scale, the 
marginally damped eddies at the damping scale will the most important. This implies a 
picture very different from the Goldreich-Sridhar cascade. 

If neutrals constitute a tiny impurity in the ionized plasma, it is clear that they cannot 
affect the MHD cascade. We shal quantify this intuitive picture below. 



2.1.3. Ion-neutral decoupling: a simple model 

Consider first a toy model of ion-neutral interaction. If we ignore viscous damping, the 
equations for the ions and neutrals are: 

— = — i^A^Vi + Fi , (12) 

'^in Pi 

and 

!!!L = I1^_1^ + F^ . (13) 

^ Pn tin 

For simplicity we have used Fi and F„ to denote both pressure forces and the nonlinear 
turbulent accelerations for the ions and neutrals, respectively. 

To linear order equations, and ignoring sound waves, (12) and (13) give the dispersion 
relation for Alfven waves in a partially ionized plasma (Mclvor, 1977): 

where is a dynamical coupling parameter, defined by 

K^^f . (15) 

Pn ''in 

When the couphng is very tight, b!; ^ 1, and we have the usual relation for Alfven and 
pseudo- Alfven modes in a single fiuid, 

T-'^-ul , (16) 



where r is imaginary. 
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As the turbulence cascades to smaller scales, r, and consequently, b!;, decreases. If 
we express the value of at the viscous damping scale as Kc, then we have two obvious 
alternatives. Either Kc > 1 and two-fluid effects are negligible right up to the damping scale, 
or < 1 and the ions and neutrals decouple in the middle of the turbulent cascade. In 
order to calculate we need to calculate viscous damping in the large K limit. 



Combining Eqs (6) and (11) we get 

Va \l \c, 



'-^-fn'(^] (k) (-]ilnk)-'/' . (17) 



where c„ is the sound speed and /„ is the neutral fraction. For most of our applications we 
will have /„ ~ 1. The damping scale, k~^, is defined by tdamp ^ Tk, so that 



.--.?/„c„.~(|)(£)"X^)l|,)^»-. (-) 

/ \ 1/2 / , N 1/2 / N 3/2 



and 



Vc ~ Vi 



CnJ V ^ / 



Prom the definition of '^c = ^("^c) we see that 

In) \PntinCnJ J \ Vi 



...... f'V'7^U-Y"7M"\ p.) 



Equation (22) seems to imply that 'i^c must always be greater than one, but a closer 
examination suggests that for pi <S p„ the third term on the right hand side can be small 
enough to offset the second term. In fact, the value of has to be determined for each 
situation. When Kg ^ 1 equations (18) through (21) have no physical meaning, since the 
plasma stops behaving as a single fluid when K drops below one. Nevertheless, we can still 
use these expressions as useful parameterizations of the turbulent cascade. In particular, if 
Kc < 1 then the scale of decoupling, when K = 1, is given by 



(23) 
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while at higher wavenumbers equation (14) becomes 

This expression has two roots. When the colhsion rate is so small that the ions and neutrals 
are completely decoupled we have usual dispersion relation for Alfven waves 

T-''^-u\p/pi . (25) 

When the collision rate is large the two roots are 

U}\tin — , - 7^. (26) 

Pi T^in 

The latter root corresponds to the case where magnetic forces are neghgible and the former is 
the usual ambipolar diffusion rate, when the magnetic field pushes the ions through a neutral 
background. Neither limit is appropriate for hydrodynamic turbulence in the neutral fluid, 
in which case r, and k are imposed by the turbulent cascade. From equation (13) we 
see that if the ions are prevented from moving by the magnetic field then we get a damping 
rate ~ {Pi/ Pn)t~in ^ which will be negligible once the hydrodynamic cascade has proceeded to 
eddies with turn over rates larger than the decoupling rate. Since the hydrodynamic eddies 
will be approximately isotropic, equation (12) guarantees that Vi <S Vn for all scales below 
the decoupling scale. 

The limit defined by equation (25) has been previously considered (sec, for example 
Kulsrud & Pearce, 1969). We will finish this subsection by considering damping in this 
regime, and the minimum neutral fraction required for neutral damping on any scale. 

It is easy to sec that equation (25) corresponds to the regime Ink ^ 1 and the waves 
are damped at a rate of t^^- Therefore, 



^'*-'H-j^fj [y;} iifKW-- (27) 

If we have a turbulent cascade in the ions, with r^i^tin > 1, then the cascade proceeds to 
very small scales, and neutral damping plays no role in its dynamics. Although the ions are 
decoupled from the neutrals, we use the usual value of Va in this expression since it appears 
only in the combination vf/{lVA), which is the energy cascade per unit mass. The ratio 
Cntin/ln IS roughly 

^ntin ^in ^i^in ^ ^in ^ Pi f (28) 

^nn Cjitni Inn Pn \ / 
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If the ions and neutral particles have the same mass, this will be ~ ^ when the neutral 
fraction is small, and a number of order unity otherwise. 

Here we have assumed that the rms velocities of the ions and neutrals are dominated 
by the thermal distribution. Our treatment should still be applicable when the turbulence 
is supersonic, but sub-Alfvenic (see discussion of the point in Cho, Lazarian & Vishniac 
2002a), as long as we replace q and c„ with the appropriate rms turbulent velocities. This is 
most likely to happen on large scales, so that /„A; ^ 1. In this case the viscosity contributed 
by neutrals stays the same so that equation (17) does not require any modifications. In the 
opposite limit, when A';/„ ^ 1, for supersonic turbulence {vfY^^ ~ V^, where VX is the Alfven 
velocity when ions and neutrals are decoupled, i.e. VX = {nnVrin + nimiY^'^ / {niUiiy^'^VA- 
However an increase in the ion velocity will result in a proportional increase in the drag 
coefficient, and a consequent drop in tin- Thus the condition provided by equation (27) does 
not change. 

Equations (17) and (27) can be used to find the minimal neutral fraction required to 
damp the turbulent cascade at scales close to Together they imply a maximum damping 
rate for kin ~ 1- We see T^Hdamp ~ 1 at that scale requires a small neutral fraction, so that 
In ~ Vntni- Assumiug a single thermal velocity we get 

/„>/„,^7xl0-(l)-'"(|)^''(^)"„-r./., (29) 

as the minimal condition for dissipating an MHD turbulent cascade through neutral particle 
viscosity. Typically, turbulent motions in the ISM are of order Cn or greater, so only a small 
neutral fraction is necessary to affect the turbulent cascade. If there is a source of noise on 
scales ^ Im such that damping given in equation (27) is ineffective, then neutral friction 
can be ignored on all smaller scales. However, we will show in the following pages that small 
scale turbulence can appear in a partially neutral plasma even in the absence of any small 
scale driving. 

2.2. The coupled regime: viscosity-damped turbulence 

In this subsection we consider strong MHD turbulence with strong neutral particle 
damping, but on scales where the one-fiuid approximation remains valid. 

In hydrodynamic turbulence viscosity sets a minimal scale for motion, with an expo- 
nential suppression of motion on smaller scales. Below the viscous cutoff the kinetic energy 
contained in a wavenumber band is dissipated at that scale, instead of being transferred to 
smaller scales. This means the end of the hydrodynamic cascade, but in MHD turbulence 
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this is not the end of magnetic structure evolution. For viscosity much larger than resistiv- 
ity, u ^ rj, there will be a broad range of scales where viscosity is important but resistivity 
is not. On these scales magnetic field structures will be created through a combination of 
large scale shear and the small scale motions generated by magnetic tension. As a result, we 
expect a power-law tail in the energy distribution, rather than an exponential cutoff. 

Here we discuss a conservative model for the damped regime. It is motivated by sim- 
ulations by Cho, Lazarian & Vishniac (2002b, henceforth CLV02b) and is consistent with 
those simulations (see also Cho, Lazarian & Vishniac, 2003b). A complete understanding of 
this regime will be deferred until higher numerical resolution runs arc available and a more 
detailed theoretical treatment has been completed. Our goal here is a model which can serve 
as an approximate guide. We will construct this model using the notion of local interactions 
in phase space (cf. Spangler, 1999), a constant cascade of energy to small scales, and a force 
balance between magnetic tension and viscous forces. We add to this two ingredients that 
are directly suggested by the simulations: a constant curvature for the field lines and an 
intermittent magnetic field distribution. 

To begin with we define a filling factor (pk, which is the fraction of the volume containing 
strong magnetic field perturbations with a scale k~^. We denote the velocity and perturbed 
magnetic field inside these subvolumes with a so that 

vl = <t>kvl (30) 

and 

hi = <t>khl. (31) 

We note that at the critical damping scale, and /cy^c; the time-scale for inviscid 
turbulence is ^ 

^ klv Ri k\\^cVA ~ KVc ~ ^^J^~yj2 ■ (^2) 

Assuming we are in the coupled regime, we have = PiTc/ {Pntin) > 1- Motions on this 
scale are marginally damped, so that smaller scale structures will be continuously sheared at 
a rate r~^. These structures will reach a dynamic equilibrium if they generate a comparable 
shear, that is 

kvk ~ T~'^- (33) 
The magnetic energy will cascade to higher wavenumbers at this same rate, so that 

blkvk ~ ^- (34) 

Consequently, bk ~ be for k > kc- 
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Next we assume that the curvature of the magnetic field fines changes slowly, if at all, 
in the cascade. This is consistent with a picture in which the cascade is driven by repeated 
shearing at the same scale. It is also consistent with the numerical work described in CLV02b, 
which yielded a constant k\\ throughout the viscously damped nonlinear cascade. 

Finally, we can balance viscous and magnetic tension forces to find 

k^uvk ~ kuT'^ ~ max[6feA;c, -Bofc||,c]^A; ^ (35) 

This suggests a picture in which the cascade to larger wavenumbers consists of an evolution 
to increasing gradients perpendicular to both the mean field direction and the local perturbed 
field component. This implies rapidly increasing magnetic pressure gradients, but these can 
compensated by plasma density fluctuations. 

Combining our results we find 



and 



h ~ he ' ~ &c, 

(kc\ fkA 

J. 



(36) 
(37) 
(38) 



The associated velocities fall rapidly, although not exponentially. In terms of one di- 
mensional spectra we have a magnetic energy spectrum E''{k) ~ k~^, and a kinetic energy 
spectrum E'"{k) ~ k~^. These scaling laws should be compared with E'" ^ ^ k~^^^ 
for the Goldreich-Sridhar predictions for the inviscid regime. The fact that the local mag- 
netic perturbations increase, albeit slowly, with increasing wavenumber implies that they 
will eventually exceed the strength of the background field. Further research should clarify 
whether this is a problem. 

If we compare this model to the simulations shown in CLV02b and the analysis in Cho, 
Lazarian & Vishniac (2003b) we see that the magnetic field Fourier power spectrum in the 
simulations is slightly steeper than expected, i.e. 

Eeik) oc k-^-^ (39) 

while the kinetic energy spectrum shows a stronger deviation 

E^{k) oc (40) 
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The fact that the deviation in the exponent is twice as large for as for Eb is consistent 
with equation (35). It remains only to explain why the magnetic power spectrum is not 
actually fiat. The most likely explanation is that since the power spectrum amplitude at a 
given k is the Fourier transform of the correlation function, it represents an integration of 
the correlation function over all scales r less than ~ k~^. Since the energy is distributed 
logarithmically on all scales between the damping scale and the dissipation scale, this in- 
troduces a factor ~ ln(/c^j„//c) to Esik). The simulations have a range ~ 20 between the 
scale of viscous damping and the scale of resistive damping, so this bends the slope of Esik) 
downward by roughly the required amount. 

The functional form of (pk is somewhat more difficult to test. In CLV03b we constructed 
filtered maps of b by separating the b(k) into broad bins. Plotting the cumulative magnetic 
energy versus volume for each filtered map gives a measure of the volume filling fraction 
as a function of scale. We show the results in figure 2 (Fig. 4 in CLVOSb). Comparing 
the concentration of magnetic energy as a function of scale wc see that, as expected, the 
concentration of magnetic power increases at smaller scales. Our results may have a slightly 
shallower dependence on wavenumber than predicted by equation (38), but this depends on 
the value of the cumulative fractional magnetic energy we choose for comparison. 

Finally, wc note that the eddies in this regime are not described by the linear solutions 
to equations (12) and (13). The transfer of magnetic power from large to small scales makes 
these solutions largely irrelevant. 

2.3. The decoupled, damped regime 

As we have seen, ion-neutral decoupling may happen either before or after the turbulent 

cascade reaches the viscous damping scale. In both cases the cascade will show similar 
behavior in the decoupled and damped regime, although the details of the onset of this 
stage will differ. We begin by discussing the case where decoupling occurs below the viscous 
damping scale, that is, after the cascade has entered the stage described in the preceding 
subsection. We will discuss the case of early decoupling, Kg < 1, in the middle of this 
subsection. We will conclude this subsection by describing the termination of the damped 
regime and the reappearance of a strong turbulent cascade involving only charged particles. 

When Kg > 1, decoupling occurs as a two step process. First, the pressure support from 
the neutral particles becomes ineffective and the energy cascade is reduced by a factor of 
{Pi/Ptot) ~ {ni/ritot)- At somewhat smaller scales the difference between the ion and neutral 
velocities becomes comparable to Vk and the neutrals can be treated as a static background. 
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We will rely on the same physical arguments used successfully in the previous subsection 
to predict this part of the spectrum. Currently numerical simulations address only larger 
scales. 

The first stage of decoupling occurs when the ambipolar diffusion rate becomes compara- 
ble to and the neutral particles infiltrate the zones of intense magnetic field perturbations. 
This happens at a wavenumber kp given by 

^in '^p ' c 

where we use the subscript 'p' to denote the pressure decoupling scale. Combining this 
criterion with equations (32) and (36) we see that 

^ ~ Hy3. (42) 

When /„ ~ 1 the streaming of neutrals across the perturbed field lines will be accompanied 
by significant dissipation, so that the turbulent energy cascade is reduced in amplitude by 
the factor by which the pressure support is reduced, that is ~ pi/p. The volume filling 
factor must also drop by this same factor in order to maintain the condition expressed in 
equation (35). In our simplified model we will model the sharp drop in bk and 0^ at /Cp as a 
discontinuity, although in reality we expect a smooth transition. It is important to note that 
there is no associated discontinuity in the rms electron density fluctuations. The volume 
averaged strength of the magnetic pressure fluctuations drops by pi/ p as the ions go from 
providing only a fraction (pi/p) of the compensating plasma pressure to supplying all of it. 
The only observational signal may be a change in higher order moments of the scintillation 
statistics. 

The neutrals and ions decouple entirely when the viscous drag coefficient becomes com- 
parable to the neutral drag term in the ion force equation. This sets in when 

pi^e ~ ^, (43) 

or at 

kd^hifn^.y/'. (44) 
For k > kd the conservation of energy condition and force balance equations become 

bl - hfj^ (45) 



and 
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This implies 



' {fnKf^ (^f J &c, (47) 

~ A^, (48) 

, \ 1/2 / , X 1/2 
"'c \ / ''in 



Equations (36), (44) and (47) imply a maximum local perturbed field strength, at wavenum- 
bers ~ kd, of 

bn.a. ~ &c(/n«c)'/'. (50) 

Under most conditions bmax will be only moderately larger than be- We note that the volume 
filling fraction increases with wavenumber in this regime. This prediction should be tested 
using a two fluid code. 

How is picture modified when we consider < 1? The decoupling condition, K ~ 1, 
evaluated to obtain equation (23), is equivalent to the pressure decoupling criterion given 
in equation (41) in this case. We see from equation (26) that turbulent motions in the ions 
are strongly damped below this limit, so that wc arc immediately in the damped, pressure 
decoupled and dynamically decoupled limit. Consequently, we expect an immediate drop in 
(pk and bk at the decoupling scale. On smaller scales we can invoke force balance and energy 
conservation to write: 

bl - b'^Pj, (51) 

and 

where r' is the eddy turn over time at the scale k' defined in equation (23). Equations (51) 
and (52) imply 

h ~ V (/,. j) ~ WK" (/»!) . (53) 

and 

where v' and b' are the rms velocities and magnetic field perturbation strengths at the scale k'. 
The differences between equations (47)- (49) and equations (53)- (55) are due to the different 
stirring rates imposed by the marginally damped eddies in the two different cases. 
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The persistence of the energy cascade to very small scales imphes that the magnetic 
field will revert to a strong turbulent cascade once colhsions with the neutral background 
are slower than the magnetic dynamical evolution rate. This sets in when 

t;:. (66) 



(47rp,)^/' 

Prom equations (47) and (53) we find that the onset of small scale turbulence is at a wavenum- 
ber 

kt-K(^j^^crmn[lX^']. (57) 

In both cases kt coincides with a filhng fraction of order unity, so that the distinction between 
bk and bk disappears. At this wavenumber small scale instabilities will cause k\\ and Vk to 
jump sharply, while bk drops. There will be a sharp mismatch between the correlation time 
of the damped turbulence on shghtly larger scales and the eddy turnover time of the strong 
small scale turbulence, which will start at ij„. This suggests that the turbulent cascade will 
be driven intermittently, with a buildup of small scale energy followed by its rapid release. 
The duty cycle for this process will be approximately the ratio of tin to the eddy turn over 
time for the marginally damped large scale or 

e~ -^min[l,«7i]. (58) 

Pn 

During the active phase of the small scale turbulence the cascade rate will be 

b{kty/Un ~ PnT^ max[l, K]. (59) 



The corresponding turbulent spectrum will be 

Vk ~ vi 



^ (kl)-'/^ (^^^ max[l, '^1% (60) 



Va 



and 



h-]m''\YA) \JJ (61) 

with hk ~ {A-KpiY/^k- 

These bursts of turbulence will terminate at some small dissipative scale, due to either 
ohmic resistivity, ion viscosity along the magnetic field lines or to plasma effects in a coUi- 
sionless medium. The resistive perpendicular wavenumber is set by the condition kres'H ~ ^k, 
or using equation (60) 



-22- 



For partially ionized gas in the interstellar medium (see Draine & Lazarian, 1998, for a 
discussion of idealized interstellar phases) the collision rate is much less than the ion cyclotron 
frequency and the ion Larmor radius is greater than k~Jg. The condition for damping due to 
parallel ion viscosity is 

klcJi > k\\VX, (63) 

where, as before, is the Alfven velocity for the ions alone, and k\\li < 1. In the limit where 
the plasma is strongly magnetized and pi <^ p„ we expect > Cn, so that this criterion is 
never satisfied. That is, damping due to parallel transport is instead given by free-streaming 
along the field lines, which has an associated damping rate which is always less than the 
cascade time scale. We have already noted that the perpendicular ion viscosity is, at most, 
only shghtly greater than the resistivity. It follows that the turbulent cascade is truncated at 
the Larmor radius, when /cr^ ~ 1. In dense plasmas, like stars and accretion disks, collisions 
dominate, ion viscosity is negligible, and the cascade ends at the resistive scale. 



2.4. Turbulence in the ISM 

Interstellar medium is turbulent with turbulence spreading over a range of scales from 
hundreds of parsec (see Lazarian & Pogosyan, 2000; Stanimirovic & Lazarian, 2001) to 
astronomical units (see Spangler 1999) As Armstrong, Rickett & Spangler (1997) pointed 
out, the fact that interstellar scintillation suggests a power law spectrum consistent with 
Kolmogorov turbulence is already a strong indication that the observed scales are connected 
by a power cascade. The lack of any feature clearly attributable to neutral particle damping 
can be seen as a counter-argument, but the model sketched above implies that it will be 
difficult to detect such a feature. For example, although the average neutral fraction in the 
diffuse ionized medium is uncertain, it probably lies somewhere between a percent and a 
few times that. From equation (29) we see that the critical neutral fraction for this gas is 
on the order of several percent, depending on the local Mach number, so fn'^c is at most of 
order unity. For low Mach numbers one might expect this to translate into an exponential 
suppression by a factor of order unity, since the eddy turnover rate will drop with the 
amplitude of motions on a given scale, allowing for more effective damping. However, since 
the electron density fiuctuations will trace the magnetic pressure fiuctuations caused by 
pscudo-Alfven modes within the turbulent cascade, the mean square variance of electron 
density on small scales, kiarge > kt, will look like an extrapolation of the same quantity from 
large scales, kgrnaii < niin[A;', fcj, but reduced by a factor of 

^ (-) (^) - {-) ^c-^^']- (64) 

\pij h[ksmallY \ksmallj \Pn J 
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This will be a very modest reduction. The scintillation power spectrum will look like a fairly 
continuous power law from large to small scales, but with a range of intermediate scales 
where it goes fiat (from kc or k' to kt) before dropping sharply to an extrapolation of its 
large scale behavior. The dynamic range of this flat region will be ~ (p„/pj) max[l, Kg]. 
Again, in the diffuse ionized medium this will correspond to a factor of order unity (or 
less) in length scales, which will not produce an observational signature. The warm neutral 
medium should show a more pronounced shoulder in the power spectrum, although given the 
heterogeneity of the local interstellar medium, detecting the signal of neutral damping using 
electron density fluctuations still represents a challenge. The coldest, and densest, phases of 
the ISM have a negligible impact on the scintillation measurements. 

In Table 1 we show critical scales and parameters for some ideahzed phases of the ISM, 
taken from Draine & Lazarian (1998). The most obvious point is that the different regimes 
discussed in this section cover a very modest range of scales. Ignoring the parallel scales, 
which may not leave a clear observational signal, the entire range of damped scales covers 
only two orders of magnitude in the cold and warm neutral phases of the ISM. For the 
bulk of the ISM, the turbulent cascade consists largely of a large scale cascade involving all 
the particles, and a smaller scale cascade involving only the ions. (We have not considered 
the role of charged grains in this analysis.) Current simulations of strongly damped MHD 
turbulence, which apply only to the range between kc and kp, already have more dynamic 
range than in the ISM itself. The small wavelength limit on ISM turbulence, set by the 
Larmor radius, is only weakly dependent on local conditions. 

In magnetically dominated environments, e.g. in molecular clouds, compressibility can 
be important. This will extend the range of scales over which the turbulence will lie in the 
viscosity damped regime. 

3. Reconnection Rates 

We arc now in a position to rederive the speed of reconnection in turbulent magnetized 
plasmas, including the effects of a neutral component. Below we briefly consider highly 
ionized plasmas, extending our previous analysis and highlighting steps in our derivation 
that will be affected by neutral particles. In §3.2 we consider stochastic reconnection in 
turbulent plasmas with a small ionized fraction. 
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3.1. Highly ionized plasmas 

For /„ < fcrit, as defined in equation (29), the ion-neutral coupling is dynamically 
insignificant on all scales. At sufiiciently small scales a modest fraction of the energy cascade 
will go into heating the neutral gas, but given the approximate nature of our discussion we 
will ignore the resulting corrections to Vk. This limit is equivalent to the fully ionized case, 
described in LV99. Here we will briefly review the major results from that paper, since 
we need to generalize them in order to calculate reconnection speeds for partially neutral 
plasmas. 

The basic geometric constraint on reconnection speeds is 

where A(L) is the width of the ejection surface for a current sheet of length L. In the absence 
of viscosity or neutral friction the ejection speed is ~ Va- The Sweet-Parker rate comes from 
taking A ~ v/^rec- In a collisionless plasma A can be ~ rx,, the ion Larmor radius, when 
this is larger than the resistive layer. Both of these estimates ignore the role that current 
sheet instabilities may play in increasing A(L). Here we are principally concerned with 
another effect, which is that as the ejected plasma moves out of the reconnection region, 
along magnetic field lines, the size of A(L) will increase as the magnetic fields lines diffuse 
away from the current sheet. 

In the presence of noise individual magnetic field lines will stay within the current sheet 
for some distance, Lmin <^ L, whose value will depend on the current sheet thickness and the 
amplitude of the noise. These individual patches of fiux will reconnect at some speed Vrec,iocai 
obtained by substituting L^j^ in place of L in equation (65). In the presence of noise, the 
current sheet contains many independent flux elements, all reconnecting simultaneously, so 
that we have an upper limit on the reconnection speed given by ~ {L / Lmin)yrec,iocai- Since 
the geometric constraint on the reconnection speed can be applied on any scale between L 
and Ljnin this implies 

Vrec < ^^^V,j,a{L'), for aU L' e [Lr^n, L], (66) 

where A(L') is the distance field lines wander perpendicular to the large scale magnetic field 
direction within a distance L' along the field lines. For a turbulent cascade the diffusion 
length A is the solution to 

^ = i..(A) (67) 

where Db defines the stochastic diffusion of field lines separated by a distance A. Assuming 
that the field line stochasticity is caused by turbulence in the medium, the functional form 
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of Db depends on the nature of the turbulence, and is hkely to be a function of scale. For 
randomly diffusing field lines this is 

Db{A) ~ max[( j kj^] for all kA > 1. (68) 

where /cy {k) is the parallel wavenumber as a function of the perpendicular wavenumber. In 
the presence of strong turbulence this expression is k~'^k\\{k) oc k"'^^^ . For long wavelength 
perturbations, kA < 1, there can be a contribution which is reduced from its value at 
the scale k~^ by factor (/cA)^. Assuming strong turbulence and taking Ak ~ 1 we see 
from equations (67) and (68) that we expect A oc k'^^^'^ ~ L'^/^. In other words, the rms 
separation of field fines should grow as the distance along the direction of the mean field to 
the 3/2 power. In figure 3 we plot rms separations for a variety of initial separations in a 
256^ cubed hyperviscosity and hyperresistivity simulation. The details of the simulation can 
be found in CLV02a. We see that the separation grows as distance along the field lines to the 
3/2 power, as expected, subject only to a time offset caused by the finite initial separation. 
Also as expected, above the dominant eddy scales the curves roll over into the usual L^/^ 
law. Independent, but similar, simulations are described in Schekochihin, Maron, Cowley & 
McWilliams (2002) who found similar results. 

In LV99, and here, we are concerned with the effect of magnetic field line diffusion 
on reconnection. However, this same effect removes the suppression of thermal conduction 
perpendicular to large scale field lines in diffuse magnetized plasmas (Narayan & Medvedev 
, 2001; Cho et al., 2003; Maron, Chandran & Blackman, 2003; Chandran & Maron, 2003). 
This has had a dramatic effect on our understanding of the thermal structure and history of 
hot plasmas in galaxy clusters, which has not yet reached a final resolution. 

For an ionized plasma with negligible viscosity, Veject ~ at all scales. Combining 
equations (6), (8), (9), (66) and (68) we find that the most restrictive limit on Vrec comes 
from L' L and gives 



Vrec < friiiin 



I ' VL 



(69) 



In a real ionized plasma the viscosity is not necessarily negligible. When the particle collision 
rate is low the perpendicular viscosity will be larger than the resistivity by a factor of several 

times the ratio of gas pressure to magnetic field pressure. Under relevant astrophysical 
circumstances this will be factor of order 10, raising the possibility that viscosity may reduce 
magnetic field line diffusion at scales slightly larger than the scale of resistive dissipation 
However, since the strongest constraint comes from the largest scales, this will not affect the 
constraint given in equation (69). 
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The claim that this is an actual estimate of the reconnection speed, rather than an upper 
limit, follows from the absence of any other important constraints on reconnection. In this 
case the most obvious alternative constraint is that individual flux elements must reconnect 
many times after their initial reconnection. If subsequent reconnection events are slow, then 
the current sheet will quickly evolve into a tangled mass of reconnected magnetic field lines. 
If the scale of the current sheet is set by the ion Larmor radius, then this will define the 
scale for all secondary reconnection events, and the small scale reconnection speed will be 
~ Va- The case where resistivity defines the width of the current sheet is more comphcated. 
In LV99 we showed that we can consider secondary reconnection events to be similar to 
the large scale reconnection, then used a self-similarity argument to show that secondary 
reconnection was fast. In either case we find that equation (69) is, as advertised, the actual 
reconnection speed for field lines in turbulent background. 



3.2. Partially ionized plasmas 

As we have shown in §2, a turbulent cascade in a partially ionized plasma is considerably 
more complicated than its counterpart in the fully ionized case. In particular, the limit 
on the reconnection speed expressed in equation (66) no longer increases monotonically 
with wavenumber, and equation (69) may not be a fair estimate of the reconnection speed. 
We need to consider the full range of dynamical regimes, and search for minima in the 
reconnection speed limit given by equation (66). This requires estimates for Veject and Db{A) 
over the full range of the turbulent cascade. We have already seen that for k < kc min[l, Kc^^] 
the most stringent limit comes from L' L, which gives a fairly generous upper limit on the 
reconnection speed. 

We begin by noting that when viscous drag is important, the ejection speed from a 
volume of thickness A and length L' is 

Veject ~ ^ max[(A;eA)^ (A^^c)"']. (70) 

At large scales, k < kc, the plasma acts as a single fluid and we have Veject ~ Ki- At very small 
scales, when the ions are completely uncoupled from the neutrals, we get Veject ~ Va{p/ Pi)^^'^ ■ 

For k > kt a, rescaled turbulent cascade emerges. As long as the current sheet thickness 
is smaller than A;^^, a rcscalcd version of our previous analysis emerges. In particular, we 
expect the most important constraint on the reconnection speed to arise from the smallest 
wavenumbers in this regime, k kt- Due to the intermittent nature of the turbulence we 
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have a factor e multiplying the usual expression. Using equations (57) and (58) we get 

n^ik > « ~ .K\.. (1)^' ^ (^) (^) ^ I (I)'' . (71) 

which is always much less than the diffusion coefficient associated with the viscous damping 
scale. However, since it applies to much smaller scales it still has physical significance. 

The diffusion coefficient for the damped scales is less obvious. The first comphcation is 
that since we are interested in the width of the ejection zone, the use of an rms value for 
A is not obviously correct. The median value would be more appropriate. Since we are not 
concerned with constants of order unity here the distinction is unimportant if the distribution 
function for field line separations is a smooth function with a single peak. Although it is not 
obvious that this condition is satisfied in the highly intermittent viscously damped regime, 
we see from figure 4 that the median and rms values of field lines separation seem to track 
one another as a function of distance along the field lines. The median is a factor of few 
lower. 

The second complication is that our model of the viscously damped regime does not 
predict the value of Dsik). Since the curvature wavenumber of the field lines is constant as 
a function of scale, we can estimate the diffusion coefficient as 

-^J (72) 

where Ic is the distance a field line remains within an intermittent structure and (pk is, as 
before, the filling factor of such structures. The value of is uncertain, but goes to k'^^^ 
for k = kc, and must fall much more steeply than k~^^^ to be consistent with the sharply 
reduced value of Dsik), relative to the undamped case, seen in figure 3. Based on a simple 
geometric picture, we will assume 

Ic^kll^. (73) 

This imphes that 

in the damped, viscously coupled regime with > 1. However, this is also the level of 
diffusion we would expect from the shear imposed at = kc- Consequently, the sharp drop 
in (pk at the scale of pressure decoupling does not affect field line diffusion. In this regime 
the value of A rises exponentially, with a constant L' ~ k'^]. We see in figure 4 that the 
simulations are consistent with comparable contributions to the diffusion coming from large 
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and small scales, and with an exponential rise in A. The sharpest constraint on Vrec will 
come from the minimum A in this regime. 

Our claim that /c^||,c -C 1 is based on the fact that the growth in field line separations 
is not consistent with the idea that a small fraction of field lines diverge sharply from their 
neighbors for a distance ~ k'^^. We can get some idea of the true situation by considering the 
full distribution of field line separations, given in figure 5. We see that while the separations 
certainly do not show a gaussian distribution, they do show that field lines join the growing 
tail of large separations at a distance which is typically shorter than k^^. 

Equations (72) and (73) imply that for k > kd the diffusion coefficient will go to a 
constant value, which persists until k — kt and (pk — 1. In this regime, kd < k < kt, we can 
use equation (57) and find 

^^(^) " • rt"'V]^ ^ ^^max[l,K,-^], (75) 
mm[A;t, k'\kc k^ 



However, this value is only relevant when, for some k > kd, it is greater than the contribution 
from the large scale shear given by ~ ^|i,cA^ (for > 1) or ~ ^\^'^ (fo^ < !)■ Once 
this value of Db applies we get A oc L'^/^ and the limit on Vrec will increase with k. The 
minimum upper limit on Vrec is at a scale 

A ^ k;' (^-^y m3x[l,K'% (76) 

with a corresponding L' of 

L' !^ k^^^max[l,^-^]. (77) 
We see from equation (70) that the corresponding ejection velocity is 

Va 

K,ect~-f niin[l,^^7i] (78) 

Jn 

Comparing equations (71) and (75) we that there is a drop in Dsik) at k — kt, implying 
that on some smaller scale where the larger scale shearing is comparable to the turbulent 
diffusion we have another local minimum in the upper limit for Vrec- This happens at a field 
line separation of 

A « k;' (^-^y (^f^y min[l,^l/'] « k;' (^-^y (^^y^'min[l,i^J, (79) 



with a corresponding L' of 



L' ^klK'-^]. (80) 
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On such small scales the ejection speed will be Va{p/ PiY^"^- 

We are now in a position to evaluate the upper limits on the reconnection speed. We 
begin with the one set at very small scales. Using equations (6), (7), (66), (79), and (80) we 
find that 

This is generally less restrictive than the limit given in equation (69). 

The limit on the reconnection speed from intermediate scales can be found from equa- 
tions (18), (20), (22), (66), (76), (77), and (78). In this case we find 

Kec < vrYjf-' min[^-, (82) 

This is potentially more restrictive than equation (69), and certainly more restrictive than 
equation (81), and needs to be evaluated for specific circumstances. This limit applies 
whenever the current sheet is narrower than the value of A given in equation (76). For the 
specific examples we consider in the next section this is usually the case. 



4. Reconnection in Veirious Phases of ISM 

All common phases of the ISM are strongly coUisionless, so that ambipolar damping 
has a reasonable chance of changing the turbulent power spectrum and consequently the 
stochastic reconnection speed. The ISM is extremely heterogeneous, but we can illustrate 
the effects of turbulence using the same idealized phases from Draine & Lazarian (1998). 
We have evaluated the reconnection speed using equations (22) and (82) for each phase. 
Our results are given in the second to last line of Table 1, and we give a detailed discussion 
below. (The last line includes the effects of tearing modes, which appear to dominate in the 
densest parts of molecular clouds). The turbulence which creates field line stochasticity is 
assumed to be supplied at some large scale, specified below, so we are implicitly ignoring 
the possibility that as reconnection proceeds it will provide a local source of turbulence^ We 
have also ignored the possibility that reconnection can change the ionization balance near 
the current sheet. 

We see from an inspection of Table 1 that neutral drag produces a limit on the reconnec- 
tion speed which is usually more restrictive than the limit given at L' — L, that is, the limit 



^The increase of stochasticity due to reconnection may be the source of the finite time instability associated 
with solar flares (Lazarian & Vishniac, 1999; Vishniac & Lazarian, 2000). 
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which applies to an ionized plasma. The only exception is the warm ionized medium, where 
the neutral content is ~ 1% and the turbulent cascade continues down to r^. Assuming 
strong turbulence, with a Mach number of order unity, the reduction in the reconnection 
speed from the ionized limit is only a factor of ~ 10 for neutral atomic gas. This is mod- 
erately sensitive to the scale of turbulent energy injection, but since we expect that L ~ / 
this will not be large effect. On the other hand, strongly subsonic turbulence implies a much 
smaller limit on Vrec- Quiescent regions in the ISM should have reconnection speeds much 
smaller than the local Alfven speed. 

Molecular gas occupies only a small fraction of the ISM's volume, but plays a special 
role as the site of star formation. Table 1 indicates that as we consider increasingly dense 
molecular gas, we move towards the case of decouphng above the viscous damping scale. 
In this limit the reconnection speed becomes more sensitive to the energy injection scale. 
We have assumed a single (large) scale here, but it may be that a much smaller value is 
appropriate in small dense regions. In that case the limit on the reconnection speed would 
drop significantly. However, if the density increases to the point where resistive instabilities 
produce current sheets substantially broader than the width given in equation (76) then 
the limit on the reconnection speed will once again increase. This accounts for the double 
entry for molecular clouds. Following the approach used in other phases we get that the 
reconnection speed is lower than the Alfven speed, even for a Mach number of 1, by more 
than an order magnitude. Unlike the atomic phases, this result scales as the square root of the 
turbulent velocity, so that weaker turbulence has a less dramatic effect on the reconnection 
speed. However, this result can be a substantial underestimate. A large scale current sheet 
is unstable to tearing modes (Furth, Killeen, & Rosenbluth (1963)), and these provide a 
localized source of turbulence which will broaden the current sheet. This sets a minimum 
value of the current sheet width. The resulting limit on Vrec takes precedence when it gives 
a larger A, and a larger Vreci than estimates based purely on the turbulent cascade. 

We show the effect of tearing modes in the last line of Table 1. For Kc < 1, for example, 
in the dense cores of molecular clouds, we can take the parallel correlation scale as 



The corresponding current sheet width, set by tearing mode instabilities is 




3/10 



(84) 



for the width of the broadened current sheet. This result is based on the tearing mode 
dispersion relation for a current sheet embedded in an neutral substrate (Zweibel, 1989), 
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and is equation (A5) in the appendix. For motion over scales comparable to /cy ^ we have 
Veject ~ Va SO cquation (66) implies 



" \tinPViJ \LVaJ \P VAtin J \Vambi J \Vambi J 

(85) 



where rjambi is the ambipolar diffusion coefficient defined by 



Vambi = —VAtin, (§6) 
Pi 



and 77 < rjamu < VaL. 



In a dense neutral medium the electron collision rate is dominated by collisions with 
neutrals, for which Draine, Roberge & Dalgarno (1983) give 

/ T 

t:^ ^ 2.62 X 10-V [j^j . (87) 

Consequently the resistivity for DC phase of the ISM is approximately 9.3 x 10*^ cm^ scc~^. 
Combining these results we obtain the last entry of Table 1. The large coefficient implies 
that the usual large scale limit, given in equation (69), will usually apply here. In the densest 
and most neutral phases of the ISM we recover the fast reconnection speeds typical of an 
ionized plasma. 

For > 1 we need to consider tearing modes in a viscous medium, with a parallel 
wavenumber ~ k\\^c- In this case we find 



Va J \Va 

which is equation (A8) in the appendix. Using this expression, and remembering that /c^i/ 
^11 cKa, we obtain 



Vree<VA{k\\,cL)i^-^J (^^J . (89) 

The resistivity for our idealized molecular cloud phase is ~ 1.3 x 10^. The resultant value 
of the reconnection speed in molecular clouds is competitive with the limit derived from 
turbulence, but is slightly smaller (meaning that the larger limit applies). Evidently tearing 
modes become rapidly more important as we go to the densest regions of molecular clouds, 
and are unimportant outside of such regions. 

Finally, we note that for coUisional fluids the magnetic Prandtl number, u/rj, is typically 
very small, so that the turbulent cascade extends down to resistive scales. In this case 
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equation (69) is our best estimate of the reconnection speed, assuming, as before, that 
Vi < Va- As an example we can consider the temperature minimum in the solar photosphere. 
Adopting T k. 4200, ~ 10^^ cm~^, and n„ w 10^^-^ cm~^, we find that v ~ c^tn ~ 
2 X lO^cm^sec"^, and 77 ~ 3 x lO^cm^sec"^. 

5. Conclusions 

Our results can be summarized as follows: 

1. In plasmas characterized by a moderate magnetic Prandtl number (resistivity less than, 
or of order, viscosity) magnetic reconnection follows the model described in LV99, so 
that reconnection is fast in the presence of noise. CoUisionless plasmas with small 
neutral fractions and collisionally dominated plasmas are typically in this class. In 
either case the turbulent cascade proceeds to either the resistive scale or the ion Larmor 
radius. For typical interstellar conditions, it is the ion Larmor radius which is larger. 

2. As the neutral content in plasma increases the turbulent cascade is interrupted at 
a scale greater than the neutral mean free path. However, the magnetic field per- 
turbations are not suppressed below this scale. On the contrary, we find that the 
magnetic field power spectrum fiattens out, and exhibits more power on small scales 
than the Goldreich-Sridhar spectrum. Velocity fluctuations in this regime are driven 
by magnetic field perturbations, whose energy is much larger than the kinetic energy 
perturbations on the same scale. The importance of this new regime of magnetic 
stochasticity goes far beyond understanding the reconnection problem. 

3. On somewhat smaller scales, once the ions and neutrals are decoupled, the turbulent 
cascade can reappear. The warm atomic gas in the ISM shows only a slight interruption 
in the turbulent cascade and the small scale scintillations are close to a continuation 
of the large scale power. However, this break will be more conspicuous in colder and 
denser phases of the ISM. In this hmit the small scale turbulence will be strongly 
intermittent. 

4. Magnetic stochasticity on scales smaller than the ion- neutral damping cut-off of the 
Goldreich-Sridhar spectrum promotes fast reconnection. The diffusion of magnetic 
field lines initially decreases rapidly below this scale, but the decoupling of ions and 
neutrals leads to a minimum value of the diffusion coefficient on small scales. The 
most stringent limit on reconnection speeds comes from considering the largest scales 
where this minimum value holds. This is in the range of the turbulent cascade where 
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ions and neutrals are decoupled, but neutral drag is still sufficient to stabilize magnetic 
perturbations. In most phases of the ISM this gives a reconncction speed which is just 
slightly smaller than the ionized plasma estimate given in LV99. 

5. Our study shows that the reconnection speed in the ISM is always much faster than the 
predictions given by the Sweet-Parker rate. Even in very dense and cold regions it is 
only slightly below the Alfven speed, and tearing mode instabilities in the current sheet 
tend to drive it back up to the Alfven speed. This process needs to be considered in 
models for magnetic flux removal during star formation, and for studies of the dynamics 
of magnetized molecular clouds. 

Finally, we note that our conclusions are moderately sensitive to the nature of field line 
diffusivity in the various damped regimes {kt < k < min[A;', fcj). We have argued that the 
shearing due to motions at /c ~ /Cc dominate until we are well into the decoupled damped 
regime. Our estimate for a constant field line diffusivity on somewhat smaller scales has 
not been tested against simulations. If this turns out to be too high (or low) then our 
reconnection speed estimate will have to be lowered (or raised) . 

AL and JC acknowledge NSF grant AST-0125544 and ETV acknowledges the support of 
the NSF grant AST-0098615. ETV would also hke to thank the hospitality of the organizers 
of the Festival de Theorie in Aix-le-Provence, 2001, where some of this work was done. 



The current sheet implicit in the Sweet-Parker picture of magnetic reconnection is well 
known to be unstable to tearing modes (Furth, Killeen, & Rosenbluth, 1963), which will 
create a turbulent zone between the volumes of unreconnected magnetic fiux. In LV99 we 
suggested that the broadening of the current sheet implied by this effect should lead to an 
enhanced reconnection rate for laminar field lines, on the order of 



Broadly similar results were obtained by Strauss (1988), except that we have calculated the 
typical size of the magnetic field fiuctuations using tearing mode theory rather than leaving 
it as a free parameter. 



A. Teeiring Modes and Current Sheet Width 




(Al) 
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Zweibel (1989) derived the dispersion relation for tearing modes in a current sheet of 
thickness A embedded in a neutral substrate. It is 



(A2) 



where 7 is the tearing mode growth rate and k\\/S. « 1. This favors modes with minimal 
i.e. transverse wavelengths comparable to the current sheet length. The resulting turbulence 
will widen the current sheet until the tearing modes are marginally stabilized by the shear 
due to the ejection of plasma, i.e. 7 ~ Va/L ^ k\\VA- 

When the ions and neutrals are decoupled the ejection velocity from the current sheet 

is 

Veject ~ ^11^ tin- (-^3) 



The dispersion relation for 7 > {pi/ p)t^^ is 



1' = (A4) 



So for the saturated modes 

This determines the width of the outflow zone if we can assume that the tearing modes 
produce stochastic mixing of the field fines, that is, that neighboring field fines are well-mixed 
within this volume. In two dimensions this is obviously not the case. In three dimensions 
it is probably valid. We are also assuming that the existence of turbulence in the current 
sheet does not lead to widespread mixing of field lines outside the current sheet. Although 
some mixing seems inevitable, the fact that the modes are driven only as long as they are 
narrow enough to fit within the current sheet, and grow only about as fast as plasma, and 
magnetic structures, arc ejected from the current sheet, suggests that any induced field line 
mixing will only spread the the outfiow zone by an additional factor of order unity. 

When the ions and neutrals are not decoupled, but viscous drag plays an important role 
we can rewrite equation (A2) as 



where we have assumed a single fiuid (7 > {pi/ p)t~^) and replaced drag by a neutral substrate 
with viscous damping. The ejection velocity in this case is obtained by balancing viscous 
drag with magnetic forces, 

- ^^eject. (A7) 
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We conclude that in this case 
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Table 1. Scales and Reconnection speeds for idealized phases of the ISM 



X 



ISM: WIM WNM CNM MC DC 

0.1 0.4 30 300 10* 

0.99 0.1 10-3 10-* 10-^ 



fc'-i 



T 8000 6000 100 20 10 

Fm 0.14 0.43 0.43 

Kc 9.l4^>l-i/3i/4 164f 7.8«3,f>!-i/3i/4 4.8/3/' A^-i/gV* Q.ZllH,^ M-"- 0^/* 

fee"' •■• 2.4 X 10"zi/4x-i/2/3-i/8 3.4 X 10i54/*>l-V2/3-i/8 8.8 X 10i34^*X-V2/3-i/8 

l^ij^^ 3.8 X lOi^Zg^^X-i^-i/* 1.0 X lOi^ig^^X-i^-i/-* 9.0 x 10i5;i/^>l-i^-i/4 

2.6 X 1013;~/'"'A4/3-i/2I 
4.0x1015/3-1/2 

9.6 X 1015z1/12;K-1/6^-5/24 1 7 X 10lB;l/12_y^-l/6^-5/24 5.2 X 10l3zl/12;K-l/6^-5/24 

fe;^! • • ■ 6.4 X 1015^-1/4 1.2 X 1015^-1/4 4.0 X 1013^-1/4 

fcfi 1.7 X 10i4>!i/2i-i/4^-3/4 13 X ioi3xi/2;,-i/4^-3/4 1.8 X 10»A4i/2i3;y*/3-3/4 2.6 X 107A4«3o''''/3-i/2 

fc-J • ■ ■ 8.5 X 10^/3-1/2 1.3 X 10"/3-i/2 1.9 x 10i3/3-i/2 4.0 x 10^/3-1/2 

ri 7.3 X 10^^1/2 1.2 X lO®/?^/^ 1.3 x 10® /^^/^ 1.3 x 10® /^^/^ 2.3 x 10® /^^/^ 

Vrec VT^^min[l,hQ/L3Q] 0.092t)TX^/3"^/'-|^ 0.077i;t>12/3-i/2 0.095i;tX^/3-i/^-|^ 0.094i;tA1-i/2/33/4-^ 

^30 ^30 ^30 ^30 ^30 

y.ec(TM) ... ... ... 0.05t;T>l-i^/®i3^"L30/3-''i/''* 21vTM-^li-^-^L3o 



k 



CO 
CO 



Note. — Here A4 is the turbulence Mach number, = vt jcn, ji = c^/Vj, -Pm is the (approximate) mass fraction contained in 
each phase, and a subscript '30' indicates that distances are given in units of 30 parsecs. The second row gives the ionization 
fra<;tion for each phase. The abbreviations in the first row denote 'Warm Ionized Medium', 'Warm Neutral Medium', 'Cold 
Neutral Medium', 'Molecular Cloud', and 'Dense Core' (in a molecular cloud). Length scales are defined in the text and are 
given here in centimeters. 
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blow up 

Fig. 1. — Upper plot: Sweet-Parker scheme of reconnection. Middle plot: illustration of 
stochastic reconnection that accounts for field line noise. Lower plot: a close-up of the 
contact region. Thick arrows depict outflows of plasma. From Lazarian & Vishniac (2000). 
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Fig. 2. — Fractional volume (X-axis) vs. fractional magnetic energy in the volume (Y- 
axis) for viscosity-damped MHD turbulence: smaller scales show a higher concentration of 
magnetic energy. 
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Fig. 3. — Root mean square separation of field lines in a simulation of inviscid MHD tur- 
bulence, as a function of distance parallel to the mean magnetic field, for a range of initial 
separations. Each curve represents 1600 line pairs. The simulation has been filtered to 
remove pseudo-Alfven modes, which introduce noise into the diffusion calculation. 
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Fig. 4. — Root mean square and median separation of field lines in simulation of viscously 
damped MHD turbulence, as a function of distance parallel to the mean magnetic field. In 
(a) we see the total result. In (b) we have removed long wavelength modes. In (c) we retain 
only the long wavelength modes. In all three cases the median tracks the mean, at a slightly 
smaller amplitude, and is consistent with exponential growth. The simulation used a 384^ 
grid. The data are based on 1600 line pairs. 
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Fig. 5. — A histogram of field line separations for viscously damped turbulence. The distri- 
bution shown here is the same one used for the curve in figure 4a for an initial separation of 
4 grid points. The tail of large separations grows smoothly and continuously, and involves a 
large fraction of all field lines despite the intermittency of the magnetic field perturbations. 



